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, , ' The semiclassical description of Skyrmions at small isospin chemical potential /i/ is carefully 

O ' analyzed. We show that when the calculation of the energy of a nucleon is performed using the 

straightforward generalization of the vacuum sector techniques (/x/ = 0) , together with the "natural" 

assumption p/ = 0{Nc), the proton and neutron masses are nonlinear in ni in the regime \fii \ < m„. 

(.~:> ' Although these nonlinearities turn out to be numerically quite small, such a result fails to strictly 

CO , agree with the very robust prediction that for those values of fij the energy excitations above the 

vacuum are linear in ni. The resolution of this paradox is achieved by studying the realization of the 

large A'^c limit of QCD in the Skyrme model at finite ^i. This is done in a simplified context devoid 

of the technical complications present in the Skyrme model but which fully displays the general 

scaling behavior with Nc- The analysis shows that the paradoxical result appears as a symptom 

^ 1' of using the semi-classical approach beyond its regime of validity and that, at a formal level, the 

D ^ standard methods for dealing with the Skyrme model are only strictly justified for states of high 

,-Q . isospin I ^ Nc- 

PACS numbers: 12.39.Dc, 25.75.Nq 
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I. INTRODUCTION 



This paper deals with a paradox associated with the description of Skyrmions [l|, y, l3| in systems with a non-zero 

isospin chemical potential /i/ 7^ — i.e., a chemical potential coupled to the third component of the isospin[^, |5|. 

f^ Before discussing this in detail it is probably of some utility to discuss why one should care about the problem. One 

IL| ' reason is simple: the problem of strong interactions with /i/ ^ is an interesting problem in theoretical physics in its 

. ^H , own right, regardless of whether it has any immediate experimental implications. From the perspective of QCD there 

j^ are certain aspects of the problem which are tractable at both small isospin chemical potentials and others which 

V;j \ are tractable at large isospin chemical potentialsjy, Q- Thus the study of this problem provides an additional avenue 

to gain insight into QCD dynamics. From the perspective of the Skyrme model — whose purpose is to mimic certain 

key aspects of QCD — the problem raises a distinct set of theoretical and mathematical issues which are interesting 

and illuminating in their own right. However, there is a more programmatic reason why this problem is of potential 

importance. One of the central problems in modern theoretical nuclear physics is the behavior of hadronic matter 

in extreme conditions of temperature or density. Ultimately, properties of such matter follows from the underlying 

theory of strong interactions, namely QCD. However, while ah initio calculations of QCD at high temperature and 

zero density are tractable via lattice simulations, the problem of dense matter is more problematic as they are afflicted 

by the notorious fermion sign problem. There has been considerable recent activity in attempting to develop viable 

lattice methods for systems at high temperature and low but nonzero chemical potential 8]. Unfortunately, these 

methods break down as the temperature decreases and are unsuitable for the study of cold dense matter. This is 
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unfortunate since this regime is of potential importance in astrophysical applications. 

Given this situation it is natural to consider models which (it is hoped) capture much of the essence of QCD, while 
being tractable in the regime of interest. One class of models which may be employed are Skyrme-type models. These 
have the virtue of encoding the scaling properties of QCD with Nc (the number of colors). Moreover, the models 
become essentially classical at large Nc and thus are far more tractable than fully quantum theories like QCD. In 
fact, over the years there has been a rather large literature dealing with cold dense Skyrmion matter [lj| 

Unfortunately, there is a fundamental difficulty in using the Skyrme model — or any other model — in place of QCD 
for cold dense matter. Even if a model is known to work well in reproducing the results of nature — and hence 
of QCD — in the vacuum sector, one is never sure whether the model is likely to continue to mimic QCD to good 
approximation in a quite different regime. Ideally one would calibrate the models against QCD in the regime of 
interest and test how they do. However, for cold dense matter this is not possible; the crux of the issue is that QCD 
is not tractable and observational data from astrophysics is, at best, both incomplete and quite indirect. 

In these circumstances it becomes exceptionally important to test models against QCD anyplace where it is both 
computationally tractable and shares features in common with the regime of interest — cold dense matter. In this 
regard QCD with /i/ =/= can play a critical role. On the one hand, unlike the case of a non- vanishing baryon chemical 
potential, the functional determinant in a path integral can be shown to be real and non-negative for systems with a 
non- vanishing isospin chemical potential Q. This means that the fermion sign problem is evaded and lattice simulations 
become practical. Preliminary simulations have already been done[lOJ. These were relatively small but increasing 
accurate lattice studies are possible 11]. On the other hand, this system has something in common with the problem 
of interest — namely a treatment of QCD with a baryon chemical potential. Thus, studies of models of QCD in the 
regime oi fij ^ can play an important role in testing models against QCD. 

Of course, if one intends ultimately to use a regime with /x/ 7^ to test a model against QCD, it is essential that 
one has computed the /z/ 7^ regime correctly in the model. It is in this context we wish to point out a paradoxical 
result that arises in the semiclassical treatment of the Skyrme model in the presence of a finite /x/ . 

The paradox concerns the behavior at small but non-zero values of /x/ and zero temperature. At zero temperature 
the behavior of QCD is known for small /i/. When j/i/l < ttt,^ the QCD vacuum state is unaltered by the presence of 
an isospin chemical potential. The only effect that adding a chemical potential has in that regime is to shift excitation 
energies above the vacuum {e.g., single particle energies for hadrons) by /ij/a where I3 is the isospin of the excited 
state. Thus in this regime all energies are strictly linear in /i/ with a slope given hy I3. At T = 0, nonlinearities only 
appear for \fj,j\ > mT^. When |/z/| reaches m^, the energy for the lowest pionic excitation reaches zero. Beyond this 
point it is energetic favorable for pions to condense. This pion condensation alters the vacuum structure and thereby 
alters excitation energies beyond the linear shift due to the direct /Z//3 contribution. 

The paradox is that calculations of the energy of a proton or neutron in a Skyrme model using the natural gener- 
alization of the techniques used in the vacuum sector yield proton and neutron masses which are nonlinear in ^/ in 
the regime /// < m,^. As we will show in Sec. Ill these nonlinearities are numerically quite small; however, the linear 
behavior in /i/ is a fundamental property of the quantum system the model is supposed to describe. Thus even small 
non-linearities appear to imply that there is something fundamentally wrong with the approach — or at the very least, 
something fundamental which is not fully understood. 

This is quite problematic on two grounds. The first is what was alluded to above. In order to use the fii ^ 
regime as a test of Skyrme type models in regimes where chemical potentials play a role, it is essential to have reliable 
calculations. The second is more vexing: if the standard approach fails at some deep level in its straightforward 
generalization from the vacuum sector to /i/ ^ 0, how certain can one be of its validity in the vacuum sector? 

In this paper we resolve this paradox — at least at a formal level. The implication of this resolution for phenomenol- 
ogy, however, remains somewhat unclear. 

The crux of the resolution involves the nature of the large Nc limit of QCD. The standard methods for dealing with 
the Skyrme model — finding a classical static soliton solution and then requantizing the collective modes to restore 
symmetries and describe physical states — are only strictly justified at large Nc. Thus in using these methods one is 
implicitly studying a large Nc world and then attempting to identify features which extrapolate back to the real world 
of Nc = 3. One feature of this large Nc world is in an infinite tower of baryon states with I — J. For Nc large but 
finite, the levels in this tower are split by effects of order l/Nc. However these effects grow with / (or equivalently 
J). For sufficiently large /, the excitation energies become too large to neglect and the semiclassical methods based 
on the 1/Nc expansion ultimately break down. The usual interpretation of the I = J tower of states for the physical 
world of A^c = 3 is that the I = J = 1/2 states correspond to the nucleon, the I = J — 3/2 states to the delta, while 
states with I ~ J > 3/2 correspond to the region beyond the scope of the approximation and are artifacts of the large 
Nc limit. 

In considering the formal large Nc limit, it is important to recall that there are critical ways in which the large Nc 
limit and the chiral limit do not commute [l3|. The standard semiclassical treatment of the Skyrme model is based on a 
pure 1/Nc expansion. Thus implicitly when studying chiral properties in Skyrme models using standard quantization 



methods one is implicitly taking the large Nc limit before the chiral limit. Thus at least formally, rriTr ^ N^, even 
though ?7T,jr is numerically small. This formal fact plays an essential role in understanding the origin of non-linearities 

for l/Lt/l < TOtt- 

The resolution to the paradox at the formal level is the following: when /i/ is treated as being of order N^ (but less 
than TO^), the natural generalization of the methods used in the vacuum is formally correct at large Nc for the state 
of the lowest free energy E — jiiI and for a band of low-lying states above it. However, these states are formally states 
of high isospin — / ~ N^ — and as such do not include the nucleon. Thus in a formal sense (based on \/Nc expansion), 
the nucleon is outside the domain of validity of the standard semiclassical methods. The non-linearities in ^i are a 
symptom of using the semiclassical approach beyond it regime of validity. 

The plan of this paper is as follows: The next section contains a brief discussion of why the hadrons' energies are 
linear in /x/ for \^j\ < rriT^. The following section contains a review of the semiclassical treatment of the nucleon 
in the Skyrme model; this yields nonlinearities in /i/ for |/x/| < to^, provided /U/ is treated as being formally of 
order N^. This is in contrast to the general results of the preceding section and as such constitutes the heart of 
the apparent paradox. In the section following, a simple toy quantum mechanical model is introduced. This model 
has central features analogous to the soliton and its semiclassical treatment but in a far more transparent form; it 
illustrates explicitly the underlying issues. The model has the virtue of being tractable quantum mechanically as 
well as semiclassically and thus enables to explicitly verify one's understanding. These model results illustrate the 
resolution of the paradox mentioned in this introduction. A final section discusses the implication of these results. 

II. THE QUANTUM MECHANICS OF A SYSTEM IN THE PRESENCE OF AN ISOSPIN CHEMICAL 

POTENTIAL 

In analyzing the Skyrme model at finite isospin chemical potential it is essential to recall that although the model 
is typically treated classically it is ultimately being used to model a quantum system — namely QCD. Accordingly it 
is essential to ensure that any treatments of the model correctly encode the underlying quantum mechanics. In this 
section, we review why elementary quantum considerations require that energy of the proton and the neutron are 
necessarily linear in /i/ below a critical value (provided that T = 0, as will be assumed throughout this paper). 

To begin let us consider the problem a bit abstractly. Let us consider a general quantum mechanical system 
with Hamiltonian, Hq- The system can be quite general with one important proviso — there must be a gap in the 
spectrum between the ground state and the lowest excited states. Suppose further that there is some operator, Q 
which corresponds to a conserved "charge" . 

[Ho,Q]=0. (1) 

For the present problem Q is the third component of the isospin (which is a conserved quantity in both QCD and in 
the Skyrme model) . Because Hq and Q commute they have simultaneous eigenstates (if there are degeneracies in the 
spectrum of Hq we will work on a basis that where Hq and Q commute they have simultaneous eigenstates) : 

HQ\n,q) ^ En^q\n,q) , Q\n,q) = q\n,q), (2) 

where q labels the charge of the state and n labels all other quantum numbers. 
Adding a chemical potential for Q to the Hamiltonian yields a new Hamiltonian 

H' = HQ-fiQQ. (3) 

This new Hamiltonian is a free energy operator; its ground state corresponds to the free energy of a system at fixed 
chemical potential and zero temperature. For field theories such as QCD (or the Skyrme model) it is the vacuum in 
the presence of the chemical potential. The commutativity of Hq with Q ensures that the eigenstates of Hq are also 
eigenstates of H' . Equation ^ implies that its eigenvalues are given by 

H'\n, q) = i?,;jn, q) with E'^^^ = (i?„,, - ^lQ). (4) 

Equation Q immediately implies that the absolute energy of any state depends linearly on fig for any value of /i/. 
It is important to note that this does not imply that excitation energies are also linear in fig for any value fig. The 
reason for this is that excitation energies are measured relative to the ground state energy and ground state energies 
need not be linear in /ig. Nonlinear behavior in the ground state energy occurs due to level crossing: the state which 
minimizes H' at one hq need not minimize it another. In the context of quantum field theories this corresponds to 
a phase transition to a new "vacuum" state with a condensate carrying the charge Q. However, provided that in a 



regime in which the ground state is unchanged from fj,Q = 0, the excitation energies are hnear in hq. For theories 
with a gap in the spectrum one expects a finite region in fj,Q over which the vacuum is unchanged and, accordingly, 
excitation energies are hnear. Now let us return to the issue of QCD with an isospin chemical potential and zero 
temperature. The proton and neutron are excitations with baryon number of unity and a third component of isospin 
of ±1/2 above the vacuum. Provided one is below the critical isospin chemical potential for a phase transition, the 
previous argument ensures the excitation energies for these must be linear in /i/. The critical issue is the value of fii 
for which the phase transition occurs. There is very strong evidence that critical isospin chemical potential occurs at 
/// = ±m^; for \fii\ charged pions condenses and the vacuum structure alters. This is precisely the picture seen in 
vacuum models based on chiral perturbation theory [6|. Thus, as a very robust prediction one expects the energy of 
the proton and neutron excitation above the vacuum to be linear in fij for \fij\ < m^. 

One might worry that there is no rigorous demonstration that the critical chemical potential is at \fij\ = 111^^. In 
order for the critical chemical potential to be smaller than this, there must be a state in QCD with mass per unit 
isospin less than m,r; this is highly implausible on its face. Moreover even if this highly implausible scenario were true, 
it would not alter the formal problem with the Skyrme model calculations. In those model calculations the vacuum 
does not undergo a phase transition — vacuum properties far from the soliton are unaltered by the chemical potential. 
Thus for the model the proton and neutron must be linear in ^/ and yet, as we will see in the following section, a 
naive extension of the standard semiclassical treatment lead to deviations from linearity. 

III. SEMICLASSICAL TREATMENT OF THE SKYRME MODEL AT /ij / 

In this section we analyze the semiclassical quantization of the skyrmion paying special attention to the issues raised 
by the presence of a finite |/i/| < 711.^. For definiteness we consider here the lagrangian of the SU{2) Skyrme model 
with quartic term stabilization and finite pion mass. It is given by 

C = -l^TriL^L-} + __Tr{[L„,L^]2} + -iJ^t:,{u + W - 2} . (5) 

In Eq.([5]), as usual, U represents the SU{2) chiral field and the Maurier-Cartan operator La is defined by La — U^daU . 
The isospin chemical potential /z/ is introduced by performing the replacement 

daU — > dalJ - i — [T:i,U] gat), (6) 

where Qap is the metric tensor in Minkowski space and ra is the third Pauli matrix. In what follows we will assume 
that 1^/1 < m^, namely that the perturbative vacuum L'^ = 1 is stable against pion condensation. 
We consider a spinning soliton configuration of the form 

U^AU{R-^r) A\ (7) 

where A and R are time-dependent isospin and spin rotations, respectively. U{r) is a static configuration. Inserting 
this form for U in Eq.fO we get 

L = -Mo + -A^fc [uja - f^iDsa) K - f^iDsb) + -Kij n, n^ + A*f K - M/^'sa) ^^■ (8) 

Here, we have used 

A-'^ A^^ UJa Ta, 

R-' R) _ ^ e^jk nk, (9) 



y 



and 



Dab = ^Tr [r. An A''] . (10) 



The static mass Mq is given by 

Mo = -y'dV[^Tr{Z,Z,} + J-^Tr{[Z„i,f} + ^^Tr{?7 + C/t-2} 



(11) 



while the tensors A^j,, A^j^ and A^j^ are the isospin, spin and mixed inertia tensors, respectively, given by 



Al 



^ii = 






cPr 



[[U,U^[Ta,U^ [u,u\n,U^] 



^Tr{Z,Z,„} - ^Tr{[L„Z,][L™,Z,]} 
^ Tr [il^p, r„]L,} + J-^ Tr { \y^[U, r,], Z,„] [U,L,r^] 






^ "^^ijk '^j 



(12) 
(13) 
(14) 



The number of independent non- vanishing components of the inertia tensors is given by the explicit form of the static 
ansatz. All the configurations to be considered below will be symmetric under 7r-rotations along any of the three 
cartesian axes. In this case one can prove that all the inertia tensors are diagonal. Thus in what follows we will 
assume that 



Ai 



A-' 



A A/ 



if a 7^ 6. 



(15) 



Now let us pay attention to the iVc-order of the different quantities appearing in Eq.®. As usual, Mq and the 
inertia tensors A's are of 0{N(,) while uja and fi^ are taken to be of 0{N~^). In order to proceed we have to determine 
which order in N^ has to be assigned to /z/. Since we are interested in values of /i/ < m^ it appears to be natural to 
take fii = 0{N^). In this case, to leading order in Nc Eq.® reads 



L(i) - -Mo + i 



a=l,2,3 



Aaa DsaDsa, 



(16) 



where Eq. (jl5p has been used. In principle the equations to determine the static soliton configuration should be 
obtained by minimizing this lagrangian. Some simplification can be performed by noting that since the chemical 
potential acts along the 3-axis in isospin space, the resulting configuration is expected to be axially symmetric along 
the third axis. In this case we have 

^33 = ^33 = ^33 = ^3, 

A(i = A^2 = A/, 
An = A^2 = Aj, 



^22 

M 



AM _ AM _ A 

I\ll — 1\22 — ^^ 



M- 



Using X)a=i 2 3 ^ZaDza = 1, it is not difhcult to find that Eq. p^ can be expressed as 

id) = -A/n 



2 A3 Ai?+ 2 (A3 -A/) ^i] SIT? p. 



(17) 



(18) 



Here, we have used the standard parametrization of the Wigner D-functions in terms of the Euler angles a,/?, 7. As 
we see, static mass Mst — — i^^^ depends on the orientation of the soliton in flavor space. However, if we assume that 
the deformation induced by the chemical potential is not too large we have that A3 — A/ << A3. In this case the 
/9-dependent term can be neglected and we have 



M., 



1 , 9 
Mo - 2 ^3 ^i- 



(19) 



Namely, at the static level we have to minimize a /i/-dependent mass which, obviously, leads to a ///-dependent soliton 
configuration. In what follows we will denote with an extra upper index the soliton quantities calculated with this 
/i/-dependent soliton configuration. Thus, the full lagrangian up to order 1/Nc reads 



L = -M, 



(m/) 



-Ay 
2 ^ 



(m/) 



{UJI - fli D31) + {UJ2 - ^J.1 D32) 



2 ./ L 1 



ni 



-A3 
2 3 



(w) 



(W3 - fli D33 + ^3) 



+A^^) [(wi - fii D31) Qi + (c^2 - M/ ^32) ^2] 



Defining the canonical conjugate momenta in the usual way, we have 

A\'"^ {iJa ~ fii D3a) + Af^ na 






J'J 



fbf 



-^ = A^^) (CO, - fii D3a) + a[I"^ n^ 



for 
for 



a 



1,2 
a = l,2 



j: 



bf 



PL 

duJa 



dL 



^3 



A^^^^ (W3 - M/ 7^33 



^3). 



(20) 

(21) 
(22) 
(23) 



Here, the upper index bf indicates that these momenta are defined in the body-fixed frame. Thus, the corresponding 

cja I'J + ^a J^J) - L reads 



hamiltonian H = Y^ , „ „ 



H 



M' 



1 

+ 2 



a=l,2,3 



D3aJ'J 



i^n 



A)> 



(^^I) 



{40'+{4'y 



^(Mi) 



rbf 



m 



2 A^^"^ 



j'/i'/ + i'/4f 



A 



(mi) 



a^a'-p') _ (Ai';^)y 



(24) 



It should be stressed that, as expected, using Eqs. pT|) - ([23|) one can easily check that to leading order in Nc the 
minimization of H leads to the same result as the minimization of L^^> . Defining T = I — J we obtain that H can be 
expressed as 



H 



M, 



(mi) 



/^/ ^3 + 



1 



A^^^^ - A 



M 



J' 



A 



ipi 



-A 



(mi) 
M 



a(mi)a(;'^) 



A 



(mi) 
M 



A 



ipi) 

M 



a(mi)a(;^) _ (a^^))' 



T^, 



P 



a: 



{fj-i 



2A 



(mi) 

M 



A 



(mi) 



A' 



(mi) 



a(mi)a(;'^) 



(AirO' 



W^ 



(25) 



where we have used the well-known relation I3 = — J2a=i 2 3 ^3aJa^ between the lab-frame components of the isospin 
and those of the spin in the body-fixed frame. It should be noted that, of course, |/ — J| < T < I + J. This means 
that for the particular case of the nucleon [J = I = 1/2) we have T = 0, 1. Since for not-too-deformed configurations 
Jy^^i , we expect that the mass of the T = 1 state will be much higher than that of the T = 0. Thus, 



it^i 



K\^'> « A 

for the physical nucleon we take T = 0. In this case, using J = I ^ 1/2, ( X 



1/4 we obtain that 



M 



(") 



M, 



it^i) -^ Ai/ 



T V + T 



2A 



{p-i 



2A 



(pi) 
M 



-Ay 



(pi) 



A 



(t^i) 



a(/^i)a'^I^i) 



A 



(26) 



So far we have considered the axial symmetric soliton configuration which is expected to be the lowest energy 
configuration in this case. As an approximation to this, in Ref. 5J a spherical symmetric hedgehog configuration was 
used. For completeness, we provide here the corresponding formulae. In such case, we have 



Ay- — Ay- — Ajj- — Ao 6ij . 
Thus, following the same steps as before we get 



(27) 



M 



C^) 



mI,^'^ t 



fj-l 



8A, 



(t^i)' 



(28) 



where Mq'' is related with AIh calculated in Ref. [5] by 



A4'''^^MH{ni) + lA[r'^ 



M?. 



(29) 



The fj,i dependence of the proton and neutron shifts (Ap and A„, respectively) of the corresponding masses from 
their values at ^7 = is shown in Fig. [T] To perform the numerical calculations we have used the standard set 
of parameters of Ref. [14l| : /tt = 54 MeV; e — 4.84; m^r = 138 MeV. We can readily see that for both the axially 
symmetry exact configuration (dashed line) and for the spherical approximate configuration (dotted line) there is 
a nonlinear behavior. In fact, the results for these two configurations are almost identical, with small differences 
appearing only for fij/rrij^ larger than 0.9. This provides an a posteriori justification for the use of the spherical 
approximate configuration made in Ref. 5]. While the deviation from linearity (full line) is numerically quite small, 
the very existence of such deviation implies an inconsistency with the quantum mechanics analysis discussed in the 
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FIG. 1: Shift of the proton and neutron masses with respect to their corresponding values at /.i/ = as a function of the isospin 
chemical potential /x/ for ^i < rrin- Full line corresponds to the result obtained under the assumption /i/ = 0{Nc^) while the 
dashed (dotted) line corresponds to the use of axial (spherical) ansatz assuming fii = 0{N^). 



previous section. The fact that the "natural" choice fij = 0{N^) leads to such an inconsistency represents the paradox 
which we would like to clarify in this article. 

One might wonder whether there is a way to perform the semiclassical quantization which leads to a linear behavior. 
In fact, this is possible provided one makes the somewhat unnatural assignment /i/ = 0{N^^). Then, to leading order 
in Nc, the lagrangian is given by 



-Mn 



(30) 



Therefore, in this case to obtain the static soliton configuration one should only minimize Mq. In this way one obtains 
a static configuration Ust which is independent of /i/ and, thus, of the spherically hedgehog type. The full expression 
for L reads 



L = -A/o + -Ao {Via +UJa- fJ-lDsa) (^a + ^a - tJ'lD^.a) 



where we note that now Ma and Aq do not depend on /i/. As before we define the momenta 

_9L dL 

^a ^ t; '-'a - 

OOJa 

Thus, the resulting Hamiltonian H = ft I — L is 



dila 



Ao(f^a 



M/ Dsa) 



H = Msol + fll Dsa Ja + 77T- 



P 

_ 

2Ao 



which leads to 
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MoT 



M/ 
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(31) 



(32) 



(33) 



(34) 



Comparing with the situation described above we see that here there is not intrinsic dependence of Mq and Ao on 
/i/. In this case we obtain the energy shift linear in ni represented by a full line in Fig[TJ Although this might be a 
resolution of the problem one is still left with the question of why the "natural" assignment /i/ — 0{N^) is not fully 
compatible with general expectation from general quantum mechanical principles while the "unnatural" assignment 
fij = 0{N~^) does. In the following section we will address this issue in detail. 



IV. A TOY MODEL 

As seen in the previous section the resuhs from the semiclassical quantization of the Skyrme model with the 
"natural" assignment /x/ = 0{N^) are inconsistent with the formal requirement of the analysis in Sect. [IT] that the 
proton and neutron energies are strictly linear in jj,j below the critical value of rn^. However, the calculation in the 
Skyrme model is technically complicated and these complications may tend to obscure the origin of the discrepancy. 
These complications are of two sorts. The first is that the system is a field theory with an infinite number of degrees of 
freedom. A few of these are collective modes which are requantized, but the remainder are internal degrees of freedom 
associated with the shape of the soliton. These internal degrees of freedom come into play when one alters the profile 
in response to the inclusion of /j,/. The second complication is that collective degrees of freedom — those which are 
requantized — do not commute with one another so the structure of the quantum mechanical system is nontrivial. 

However, the general scaling behavior of the system with Nc — which we believe to be the core of the problem — can 
be illustrated in a much simpler context: a quantum system with two degrees of freedom; one which plays the role 
of the internal degrees of freedom and the other the collective ones. This system is supposed to be analogous to 
the soliton excitations acting above the vacuum and hence the analogy only holds in the regime where the isospin 
chemical potential has not induced a phase transition of the vacuum — which is the regime where our paradox has 
been identified. The absolute energies in this model thus correspond to the excitation energy of the baryons above 
vacuum. The collective degree of freedom in the model is associated with a conserved quantity (which we will take to 
be the analog of the third component of isospin); we can add a chemical potential for this quantity. The system can 
be chosen so that there is a direct analog of Nc which controls the extent to which the quantum mechanical system 
can be treated semiclassically. Such a "toy model" is a useful place to test approximation schemes since it is solvable 
at the quantum level and, thus, one can directly compare approximate with exact results. 

The model, in units with fi = 1, is given by the following Hamiltonian which corresponds to a particle of unit charge 
moving in two dimensions: 
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p~A{r) +\V{\-^r'^) (35) 



where A plays the role of Nc and A is an infinitely thin magnetic flux localized at the origin and with a strength of 
precisely one half of a fiux quantum: A(x) = — ^^^^^^ where x and y indicate unit vectors (as opposed to quantum 
operators). The purpose for including this vector potential will become clear presently. Note that this model is axially 
symmetric and hence commutes with L, the two-dimensional angular momentum; the Hamiltonian can be written as 

H = Hr + 2mr'^ whcrc Hr acts entirely on the radial degrees of freedom and the L — 1/2 structure refiects the 
presence of the vector potential. It is convenient to introduce a new operator I = L — 1/2. Note that a well-defined 
Hilbert space {i.e., a single-valued wave function) requires that the eigenvalues of L be integers; accordingly the 
eigenvalues of / are ±1/2, ±3/2, ±5/2 • • •. Thus / plays the role of I^ in the Skyrme model and the vector potential 
was included to ensure half-integral values for I. 

Note that the model has a magnetic flux localized at the origin. Accordingly, the vector potential is undefined at 
the origin. Thus, to be sensible our model must exclude the particle from hitting the origin. This can be achieved by 
choosing a form for V which is repulsive at the origin and sufficiently singular. Having a singular repulsive potential 
at the origin has another virtue: if there is attraction elsewhere in the potential, there will be a minimum of the 
potential as a function of r. A classical solution localized (in r) at this minimum plays the role of the soliton; angular 
motion is then the collective motion. 

A convenient choice for the functional form of the potential is V{r'^) — ^^^.^ + ^rauP'r^ , since this system is exactly 
solvable. Including our scaling factor A (playing the role of Nc) the effective potential for radial motion associated 
with a state of fixed / becomes 



and the radial Schrodinger equation becomes: 



11^ X M+X rnuj'^i 

-7. drrdr + A \ + —— 

Zm r \ Zrav^ l\ 



^i{r) = Eii;{r) (37) 



This can be solved exactly. The energies of the lowest-lying state for any / are given by 

Ei = uj(l + v/A2 + p) 

= ^(a + 1 + ^ + 0(A-2)^ . (38) 

Note that these energies at large A accord with our large Nc expectations: the overall energy of the states go as 
A '^ A'c as do the mass of baryons while the energy splitting between states with different P (of order 1) scales as 
A^i ~ N~^ as does, for example, the N-A mass splitting. 

Before putting a chemical potential on this system, it is useful to show that a semiclassical approach analogous 
to that used from Skyrmions reproduces the leading large A (large Nc results) for both the overall energy and the 
splittings. As a first step one needs to find the analog of the soliton: a static solution of the classical equations of 
motion. This amounts to finding rg, the value of r which minimizes the potential. It is straightforward to find tq and 
the value of the potential at its minimum, Eq 



A 

V muj 
Eq = Aw (39) 

Note that Eq is the analog of the "static soliton mass" Mq appearing in e.g. Eq.®; a direct comparison with 
Eq. ([38]) shows that it accurately predicts the leading contribution to the energy. There is collective angular motion 
around the minimum. The corresponding moment of inertia is Aq = A(ro) with A(r) = mr^. It is easy to see that 
the requantization of the collective motion — in analogy to the standard treatment of the Skyrme model — yields a 
collective motion contribution to the energy of 






ET" - 77^ = TT^ (40) 

^ 2mr2 2A ^ ' 



which directly reproduces the leading /-dependent contribution of Eq. (jSSp . Thus the analog of the standard Skyrmion 
semiclassical approach correctly describes the dynamics of this toy model up to the accuracy for which it is supposed 
to work. 

In doing this comparison to Eq. (|38p the classical "soliton" correctly gives the order \ ^ Nc part of the exact 
energy. However, there is also a term of order unity. In the Skyrme model terms of order N^ correspond to the zero 
point fluctuations of the pion fields about the soliton background. In Skyrme models they are difhcult to compute and 
generally neglected. By analogy, in this model the A^ piece of the energy should correspond to the zero point motion 
of the internal degree of freedom. This zero-point energy is given by 1/2 the small amplitude vibration frequency 
about the minimum: 

1 / d-^V 

-- uj . (41) 

It is worth noting that this exactly reproduces the term of order A° ^ iV° in Eq. ([55)1 

Now consider what happens if we add a chemical potential for / to the Hamiltonian: H' = H — jiiI. The exact 
eigenstates of H' corresponding to the lightest state with fixed / is given by Ej{fii) = Ei — /x/J. Since we are 
interested in the effect of the chemical potential, it is useful to focus on the difference of the energy from the ground 
state "baryon" at /i/ = 




AEii^ii) = E'ji^ii) - Ey^ = cj v/A2 + /2 - Wa2 + - - fiji , (42) 




which is plotted in Fig. ^ for a number of different values of / for the case A = 15. 

We can now proceed to treat the problem semiclassically using an algorithm analogous to the one in the previous 
section. If one implements this assuming /i/ ^ 0{l), one has a "soliton" whose static energy, in analogy to Eg. p^ . 
is obtained by minimizing 

V.ir,,,) = Vlrnr)4 Mr) ,j ^ A ( -^ + !^) - '-mr^ ,j . (43) 



.cj, V ,~-^ .H' >^^i~-\2^^ -^^ - 2 
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{a)AEj versus chemical potential 



a 3 




(b)The small fij regime 



FIG. 2: Energy shifts as defined in Eq. (|42|l (in units of uj) versus tlie chemical potential /ij for A — 15. The solid lines represent 
exact solutions. At fij — the lowest line corresponds to 7 = 1/2, the line above it to / = 3/2 with / increasing by one unit 
each line up. The dashed line corresponds to the semiclassical solution of Eq. (|44|l Figure (a) goes out to ni — .6 and illustrates 
how the semiclassical result follows the minimum of the exact solution. Figure (b) shows the low ni region and illustrates the 
breakdown of the semiclassical treatment at very small ni . 



Minimizing this with respect to r to obtain the analog of the sohton hegdehog mass Mnipi) (see Ea. (|29)) ) one gets 



>i) 



Ai/2 



ml/2(,^2„ ^2)1/4 



EH{^iI) - K.t(r^\/i/) = A^c.2_;,2 



(44) 



where /S.Eh{ij.i) is the shift in the classical "hedgehog mass" due to the chemical potential. It is also plotted in 
Fig. ^ for the case A = 15. 

It is apparent from these plots that the classical treatment of the "soliton" works in the sense that it describes the 
shift in the minimum energy solution as a function of /i/. Moreover, it is equally clear how it works. Although each 
quantum level is linear in /i/, as /i/ increases the quantum levels cross and the value of / of the lowest state increases. 
The classical calculation tracks this level crossing — in a continuous way. Thus it is able to accurately reflect the 
behavior of the minimum energy state when /i/ is large enough so that / 3> 1. Parametrically, this regime occurs for 
^i /uj ^ 1/A ^ 1/iVc where, therefore, rotational corrections can be neglected. Conversely at small /i/ one expects 
the semiclassical analysis to fail to reproduce the shifts — and, as it can be seen in Fig. 2b, this is indeed the case. 

Let us now look at the analog of the semiclassical calculation of the proton's energy. Following the same steps as 
in previous section we get for the present toy model 



E' 



E, 



{pi 



2 



1 



(45) 
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FIG. 3: Energy shifts for the "proton" (/ = 1/2) state as a function of /i/ for the case of A 
of u). The sohd hne is the exact result; the dashed Une is the semiclassical resuh. 



15. All energies are given in units 



where 






.J^il^i 



M? 



A(''^) = m r^^^^ 



The explicit form of "proton" energy resulting from this calculation is given by 



E;'{^^I 



A"! 

2^^ 



1- 
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M? 



2u; 



(46) 



(47) 



For ease of comparison we introduce the shift in the semiclassically calculated energy from fij = defined by A!f — 
Ep'^ipi) — Ep'^{0). The corresponding exact result from Eq. ([i^ is A'^^^'^^ — ~fJ-i/2. Both the exact and semiclassical 
result are plotted in Fig. [3l It is apparent that the semiclassical procedure analogous to that in the preceding section 
does not accurately reproduce the exact result except at very small /i/ where nonlinearities are negligible. The disease 
it suffers from is precisely the one identified for the analogous Skyrmion calculation. Namely, that in contrast to 
the exact result, the calculated energy is nonlinear in fij with a nonzero curvature all the way down to fij = 0. It 
is also clear why it fails. The "soliton" used in the calculation was the minimum classical solution with fixed ^/. 
Now it is clear from Fig. [21 that for /i/ of order one (A°) this classical configuration corresponds to quantum states 
with relatively high /, i.e., states which parametrically have / ~ A. In situations where a classical solution breaks 
a symmetry, and the underlying quantum system is in the semiclassical limit, a classical configuration is associated 
with a band of states which are all qualitatively similar and all have energies which are near this minimum. Thus, 
when /i/ is parametrically of order unity, the classical solution is valid only for describing states whose energy differs 
from the minimum by energies of order 1/A. The "proton" state is not in this class. 

To summarize the results of a semiclassical treatment of this toy model for /z/ ^ 0: in the large A limit (analogous 
to the large Nc limit) and fu of order unity, the analog of the classical soliton gives an accurate description of the 
lowest-lying "baryon" states in terms of free energy — states which have / ~ A — but does not accurately describe the 
"nucleon" state. 



V. DISCUSSION 



The toy model of the previous section illustrates the resolution of the paradox of the existence of small nonlinearities 
in the nucleon's energy as a function of /i/ in Skyrme model calculations. The key point is that at large Nc and zero 
temperature with fij r^ N^ but below the phase transition (i.e., \fj,j\ < m.,^), semiclassical calculations are formally 
valid for the lowest-lying baryons states (which have / ^ Nc) but not for the nucleon. 

This formal resolution of the paradox is useful for understanding what is happening from a mathematical perspective. 
However, important phenomenological problems remain. The crux of the issue concerns the fact that in nature Nc = 3 
and the parameters of the model are fit from the Nc = 3 world. The high isospin states critical to the formal resolution 
of this problem do not exist as baryon resonances in the physical world. Thus, there appears to be no useful regime for 
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which the semiclassical method for describing how baryon free energies vary with fij . The issue is further comphcated 
by the fact that m^, while formaUy of order N^, is, in fact, very smaU due to approximate chiral symmetry. Indeed, 
rriTr is smaller than Ma — M^ which is formally of order l/Nc- Thus if one were to increase /z/ starting from zero, 
even before the first level crossing occurs in the baryons (the I3 — 3/2 state of the A crossing the proton), a phase 
transition occurs and the semiclassical methods discussed here need to be modified. It is not immediately obvious 
how to do this. 

Apart from phenomenology, an interesting formal question remains: If the semiclassical method fails, how does one 
accurately describe the nucleon (and other low isospin states, such as the A) for /// 7^ in the large Nc limit of the 
Skyrme model)? In the regime |/zj| < m^, where the paradox discussed here is manifest, this is easily accomplished. 
One first does semiclassical quantization at /^/ = and identifies the quantum states. Subsequently one imposes the 
chemical potential at the quantum level on these states. Such a prescription is guaranteed to give the correct linear 
behavior. However, the problem gets more interesting when |/i/| > ttt.^. In this regime pions condense. It is clear that 
new methods need to be developed to describe baryons in this regime. 
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